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The decay rate 7 of an excited dipole molecule inside a waveguide is evaluated for the strongly 
coupled matter-field case near a cutoff frequency u c without using perturbation analysis. Due to 
the singularity in the density of photon states at the cutoff frequency, we find that 7 depends 
non-analytically on the coupling constant g as g 4 ^ 3 . In contrast to the ordinary evaluation of 7 
which relies on the Fermi golden rule (itself based on perturbation analysis), 7 has an upper bound 
and does not diverge at ui c even if we assume perfect conductance in the waveguide walls. As a 
result, again in contrast to the statement found in the literature, the speed of emitted light from 
the molecule does not vanish at uj c and is proportional to cg 2 ^ 3 which is on the order of 10 3 ~ 10 4 
m/s for typical dipole molecules. 
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It has been well known for many years that enhanced 
and inhibited spontaneous emission can be observed for 
an atom in a resonator above and below a cutoff fre- 
quency 0,0- An explanation of this alternation has been 
proposed by Kleppner based on the combination of the 
Fermi golden rule for the transition rate and a change of 
the density of photon states in a resonator. In his calcu- 
lation, the singularity in the density of states at a cutoff 
frequency uj c may lead to a vast increase of the decay rate 
7 of an excited atom or molecule compared to the free- 
space value. This prediction has been verified by several 
experiments 0-0 ■ However, Kleppner's approach, while 
going in the right direction, is inconsistent. This is ap- 
parent because the Fermi golden rule is not applicable in 
the domain of a vast increase in the decay rate, as this 
rule is only valid for a weakly coupled matter-field system 
where the perturbation analysis may be applied. Due to 
this inconsistency, 7 diverges in Kleppner's theory for the 
case in which the resonator is made of a perfectly con- 
ducting material and the characteristic frequency oj\ of 
the atom approaches uj c - The decay rate 7 then vanishes 
for lo\ < uj c . 

The purpose of this Letter is to re-evaluate 7 for a 
strongly coupled case around oj c without using perturba- 
tion analysis. The result shows that 7 is finite at uj c even 
if we assume perfect conductance. We find that the maxi- 
mum value of 7 depends non-analytically on a dimension- 
less coupling constant g as g 4 ^ 3 (instead of g 2 which is 
assumed in Fermi's golden rule). For small g <C 1, 7 is 
enhanced by a factor of g~ 2 ^ 3 over the free-space value 
of 7 that is proportional to g 2 . Moreover, the frequency 
value at which the unstable solution disappears is found 
not at lui = lu c , but at uj\ = uj s , which lies slightly below 
lu c with uj c — uj s ~ g 4 / 3 uj c . Again, this value is much 



larger than the value ~ g 2 predicted by the perturbation 
analysis for the system where the density of states has 
no singurality at io c . 

To present the results, let us consider a dipole molecule 
(such as HC1, NaCl or KBr) that is on the order of a 
nanometer in size and which has a charge of +Ze on 
one end and charge — Ze on the other end. This charge 
couples the dipole to the electromagnetic field in the rect- 
angular waveguide. The waveguide runs parallel to the 
z-axis and extends to infinity in cither diection. We refer 
to the width of the waveguide in the x-direction as a and 
to the height of the waveguide in the y-direction as b. We 
assume b > a. The origin of z is chosen at an arbitrary 
point along the infinite waveguide. We choose the origin 
for x and y at the lower-left corner of the cross-section of 
the waveguide. We assume that the center of mass of the 
dipole is initially located near the center of the waveguide 
at coordinates (a/2, 6/2, 0). The characteristic frequency 
of vibration ui\ for a typical dipole is on the order of 10 14 
to 10 13 Hz, corresponding to a wavelength A with 10 to 
lOOyum, which is in the infrared range. Since the width of 
the waveguide is much larger than the size of the dipole, 
we can neglect the forces acting on the dipole from the 
walls of the waveguide. 

Under these conditions, one may write the Hamiltonian 
of the system with the reduced mass fi = m a mb/(m a + 
nib) for the dipole and with Z a = — Zb = Z as 
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where r = (x,y,z). The vector potential inside the 
waveguide consists of the TE modes (Transverse Elec- 
tric field modes) and TM modes (Transverse Magnetic 
field modes), A = A TE + A TM S\. For the case of 
the waveguide made of a perfectly conducting material, 
the solution of the sourceless Maxwell equations with the 
Coulomb gauge inside the rectangular waveguide leads to 
these modes in terms of the normal coordinates of the 
fields q a k (with a = E or M and k = (to, n, k) with to, 
n integers and the continuous variable k) as 
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where Wi t k(r) = cos(mirx/a) sm(mry /b) exp[ikz], 
W2,k(r) = sin(TO7ra/a) cos(mry/b) exp[ifcz], W^3,k(r) = 
sin(TO7rx/a) sm(niry/b) exp[ikz], and Bi are unit vec- 
tors. We put C mn = 2c 3 / \aba 2 nn ) with a mn = 
^(mir/a) 2 + (nn/b) 2 , uj^ = c^Jk 2 + a 2 ln , J s dk = 

Sm,n>0 I^oo and Imn = 1 for 771 and 71 ^ °J = 

1/2 for m or n = 0. 

We may rewrite the Hamiltonian in terms of the rela- 
tive coordinate ri = r& — r a , the center-of-mass coordi- 
nate R, and their canonical conjugate momenta p! and 
P. In this discussion we are interested in distances much 
larger than the size of the molecule, so that the particle 
interacts with the field at approximately the center-of- 
mass (the so-called dipole approximation). We assume 
that the velocity of the center-of-mass is so slow that 
we can neglect its kinetic energy. We also assume the 
field is weak enough that we can neglect terms that are 
proportional to A 2 . We first consider the case in which 
the dipole oscillates in the x direction (thus i"i = x\e x 
and pi = p\e x ). The extension to arbitrary direction 
will be discussed later. We then introduce the unper- 
turbed normal coordinate qi = fuo\/2\x\ + i(p±/ /jwi)] 
of the dipole. With these assumptions, the Hamiltonian 
|JT)| may be approximately written in a bilinear form with 
respect to the annihilation and creation operators (which 
are related to the normal coordinate through q a = ^/ha a ) 
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with the dimensionless coupling constant g = 
^(Ze^ui/Onc 3 ). We have g ~ 1CT 6 to 1CT 7 



for the typical dipole molecules. The interactions 
are given by Ve,)s. = -i{n/b)y/ I mn F ltk (R), and 
Vmm = -(md/au) k )Fi, k (R), with i*\,k(R) = hir 2 C mn 
Wi t k(B.)/ y/iDb.. The operators a a satisfy the usual com- 
mutation relations. For each (to, n) mode, the continu- 
ous spectrum Wk is bounded from below at ca mn . These 
lower bounds form a set of cutoff frequencies in the 
sense that only electromagnetic modes with frequency 
t^k > ca mn may propagate inside the waveguide. Among 
these branches of continuum, the TE mode with to = 
and n = 1 has the smallest value for its cutoff frequency 
u> c = caoi = 7rc/6. 

The Hamiltonian 10} has exactly the same structure 
as the Hamiltonian for the well-known Friedrichs model 
with virtual processes [ToL Hi) . This model has been 
investigated extensively in order to analyze the spon- 
taneous decay of an excited atom ^(J~Ei|- Since the 
Hamiltonian is bilinear, one can find its exact diagonal 
form as H = Tiuixb^bx + J s dk huj^b^b^, where uj\ is 
the shifted real frequency for the stable dressed harmonic 
oscillator. The new dressed annihilation operators b a are 
obtained by the Bogoliubov transformation and satisfy 
the usual commutation relations. The explicit form of 
the Bogoliubov transformation is presented in 0, 0] . 

In this short letter we limit our studies to the shifted 
frequency and the decay rate of the harmonic oscillator. 
These are found by solving the dispersion equation £ (z) = 
for the dressed harmonic oscillator, where 
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The solution of £(z) = gives a real solution z = zq = 
toi with a shifted frequency for the stable mode of the 
harmonic oscillator and a complex solution z = Z\ = 
u>i — vy with shifted frequency lj\ and decay rate 7 for 
the unstable mode. 

To find the explicit form of these solutions, let us con- 
sider the case in which ui\ is located below the next small- 
est cutoff frequency u> C 2 . The geometry of the waveguide 
can be chosen in such a way that the cutoff frequency 
w c is well separated from lu C 2 , and as such, we may have 
^1 < w c , or uj c < u>\ < uj C 2 where u>x is much closer to 
uj c than lu c2 . For this case the predominant contribution 
to the integration in Eq. (J5J comes from the component 
with a — E, to = and n — 1. Approximating the inte- 
gration by leaving out all other components and explic- 
itly performing the integration, we obtain the dispersion 
equation 
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where we have introduced the dimensionless variables, 
w\ = loi/uJc, C, = (z/loc) 2 , and a dimensionless constant 
G 2 = (4:c/auj c )sm 2 (TTRy/b). The minus branch of Eq. 
© gives the real solution z — zq for the stable mode, 
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while the plus branch gives the complex solution z = Z\ 
for the unstable mode. In the following discussion we 
shall consider the case where the dipole remains near the 
center of the waveguide, without loss of generality, and 
put Go = G\ji —b/2- Squaring Eq. |JB}, we obtain a cu- 
bic equation for £. Using the standard method to solve 
the cubic equation, one can explicitly find the shifted fre- 
quency and the decay rate. However, since the equation 
was squared, we must take care to exclude non-physical 
solutions. The discriminant of the cubic equation is given 
by D( Wl ) =4p 3 Wl 



1w l > where 
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The stability of the dipole is determined by the critical 
frequency at w\ = w s , which is given by D(w s ) = 0. This 
gives 



= 1 - 
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For g < 1 we have uj s ps uj c [1 - 2 1 / 3 (3/4)(gG ) 4/3 ] where 
lu s = w s oj c . For a given value of lu c only the stable mode 
exists with real z = Zq for tui < lu s , while both the stable 
mode and the unstable mode with a complex z = Z\ exist 
for loi > uj s . 

Then, we define a w ± = (~q w ± \J + 4p 3 )/2 in order 
to write the solutions of Eq. © as 
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where n — for the stable solution and n = 1 for the 
unstable solution. 

In Fig. 1 we plot the value of w\ — u)i/u> c and w\ = 
u!i/u> c — He(zi/u> c ) as functions of wi = u>i/ui c with a 
fixed value of u> c . In this and the next figures we set 
Go = 2 and use a large coupling constant g = 0.005 to 
exaggerate the effect of the interaction. We indicate the 
location of the critical value w s . The thin line is y = w\ . 
The thick curve below the thin line is wi for the stable 
mode. The thick curve above the thin line is w\ for the 
unstable mode. The two dashed curves in the domain 
wi < w s are unphysical real solutions obtained from z n 
by putting n — 1 and 2. 

We note that the stable solution w\ exists for all values 
of wi inside the waveguide. This is the result of the 
singularity in the density of states at the cutoff frequency 
u>k = u c in Eq. JSJ appearing as cdk — uj^duik/ {wf, — 
uil) 1 / 2 . Because of this singularity, there may be a large 
deviation in the value of w\ from w\ for ioi > 1 no 
matter how small the coupling constant g might be. This 
is a striking difference from the ordinary Friedrichs model 
which is used to analyze the spontaneous emission of the 
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FIG. 1: The vertical axis represents the shifted frequencies 
wi = loi/loc and W\ = u>i/w c , and the transverse axis is wi = 
uji/uj c with the value of ui c fixed. We indicate the location of 
the critical value w s - The thin line is y = Wi. The thick curve 
below the thin line is wi for the stable mode. The thick curve 
above the thin line is wi for the unstable mode. The two 
dashed curves are unphysical real solutions which correspond 
to n — 1 and n — 2 in Eq. © for W\ < w s . 



photon from the atom located in the vacuum without 
boundary [l7| . In contrast to w\, the deviation of w\ 
from w\ for the unstable mode is always small for g <C 1. 

In Fig. 2 we plot the decay rate j/w c — Im(zi/ti; c ) of 
the unstable mode for a fixed value of uj c as a function of 
wi. The maximum value of 7 in our system is obtaind 
at Wi = 1 (i.e., uj\ = uj c ). At this point we have p Wl — 
(Xw!+ = 0, q Wl = a Wl - = (<?Go) 4 , and we obtain for 
<7< 1, 
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The maximum value is a nonanalytic function at g = 0, 
and hence one cannot obtain this result from perturba- 
tion analysis. For j « 1 this is much learger than the 
decay rate 72 = 2g 2 (c/a)(u> c /uji) which is found in the 
perturbation region where |wi — l|/|w s — 1| 3> 1. Indeed, 
the enhanced factor for the decay rate "/ m ax/l2 oc g~ 2 ^ 3 
is extremely large for the case wi > 1 with wi ~ 1 and 
9< 1. 

Notice that the critical value w s of the unstable mode 
is located below the cutoff frequency, i.e., w s < 1. This is 
another striking difference from the ordinary Friedrichs 
model, in which the critical value is located inside the 
continuous spectrum of the field 01 • 

We note that 7 does not depend of the size of the 
molecule in the dipole approximation. One can show 
that the vibrating motion of the dipole is stable when it 
oscillates in the y direction, or in the z direction for our 
illustrated case with cjj < ui C 2- As a result, we have the 
same value of 7 for molecules oriented in any direction of 
space. 
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FIG. 2: Decay rate 7/w c of the unstable mode as a function 
of wi — u)\/u c for a fixed value of uj c - The decay rate does 
not diverge at wi — 1 even though we have assumed the 
waveguide is made of a perfectly conducting material. At 
lifi = 1 we have the maximum value of 7, which depends non- 
analytically on the coupling constant g as presented at Eq. 

E3). 



We also note that the g 4 ^ 3 law obtained in our sys- 
tem is rather universal around the cutoff frequency, as it 
can be shown to be independent of the particular model 
we choose. Indeed, one can find the g 4 ! 3 law by esti- 
mating the lower bound of u)\ at which the Fermi golden 
rule is applicable in the vicinity of the singularity in the 
density of states at u c . Due to lack of space in this Let- 
ter, we will present this estimation elsewhere. Finally, 
we remark on the velocity Vf of the light emitted from 
the dipole. Through the exact form of the Bogoliubov 
transformation, we find 




In the perturbation region \wi — l|/|u; s — 1| 1 in which 
|cDi — ~ g 2 coi and 7 ~ g 2 u>i , this expression reduces to 
the well-known formula for the ordinary group velocity 
Vf ~ v g = c[l — (cJc/cji) 2 ] 1 / 2 inside the waveguide [8j. 
However, since the perturbation analysis fails at io\ = 
lu c , one should not conclude from the form of v g that 
the group velocity vanishes at the cutoff frequency as 
stated, for example, in 8]. Indeed, Eq. lfTT|) leads to 
Vf « 2cg~ 2 / 3 ^f max w (v / 3/2)g 2 / 3 c at wi = lo c for the 
group velocity. For typical dipole molecules, Vf at uj\ = 
oj c is on the order of 10 3 ~ 10 4 m/s, which is comparable 
with the speed of phonons propagating on the walls of 
the waveguide. 
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